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Abstract 

A study of the d-dimensional classical Heisenberg ferromagnetic model in the presence of a 
magnetic field is performed within the two-time Green function's framework in classical statistical 
physics. We extend the well known quantum Callen method to derive analytically a new formula 
for magnetization. Although this formula is valid for any dimensionality, we focus on one- and 
three- dimensional models and compare the predictions with those arising from a different ex- 
pression suggested many years ago in the context of the classical spectral density method. Both 
frameworks give results in good agreement with the exact numerical transfer-matrix data for the 
one-dimensional case and with the exact high-temperature-series results for the three-dimensional 
one. In particular, for the ferromagnetic chain, the zero-field susceptibility results are found to 
be consistent with the exact analytical ones obtained by M.E. Fisher. However, the formula de- 
rived in the present paper provides more accurate predictions in a wide range of temperatures of 
experimental and numerical interest. 

PACS numbers: 75.10.Hk, 05.20.-y, 75.40. Cx 
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I. INTRODUCTION 



Classical spin models are widely studied in statistical mechanics and play an impor- 
tant role in condensed matter physics but also in other disciplines, such as biology, neural 
networks, etc. Despite their intrinsic simplicity with respect to the quantum counterpart, 
classical spin models show highly nontrivial features as, for instance, a rich phase diagram 
and finite-temperature criticality. 

Remarkably, in many relevant situations, the investigation of classical spin models allows 
to obtain several information for realistic magnetic materials. Indeed, they have turned 
out to be extremely versatile for describing a variety of relevant phenomena. This justifies 
the significant effort which has been devoted for optimized implementations of Monte Carlo 
simulations of spin models. Besides, there are several questions related to classical spin 
models which may further stimulate in developing new available computational resources, 
more efficient algorithms and powerful techniques for obtaining satisfactory answers. 

Along this direction, several methods have been employed to investigate different classical 
spin models such as Ising, Potts models and several variants of the basic Heisenberg model. 

A very efficient tool of investigation, in strict analogy of the quantum many-body tech- 
niques, is constituted by the two-time Green functions (GF) framework in classical statistical 
physics, as suggested by Bogoliubov Jr. and Sadovnikov [1] some decades ago and further 
developed and tested in Refs. [2-5]. 

The central problem in applying this method to quantum jo] and classical {2-5, 7| spin 
models is to find a suitable expression for magnetization in terms of the basic two-time 
Green function for the specific spin Hamiltonian under study. 

The same problem occurs in the related quantum {sl and classical [l-5| spectral density 
methods (QSDM and CSDM) which have been less used in literature although they appear 
very effective to describe the macroscopic properties of a variety of many body systems j^. 

For quantum spin-1/2 Heisenberg models, exact spin operator relations allow to solve 
directly the problem. The case of spin-S was solved in an elegant way by Callen [9j (in 
the context of the two-time GF method) providing a general formula for magnetization, 
successfully used in many theoretical studies. 

Unfortunately, for classical Heisenberg models, a similar formula is lacking in the classical 
two-time GF framework [l-5|. This difficulty is related to the absence of a kinematic sum 
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rule for the 2;-component of the spin vector as it happens in quantum counterpart. 

Almost three decades ago in the context of the CSDM , a suitable formula 
for magnetization was suggested for the classical isotropic Heisenberg model. However, its 
analytical structure was conjectured on the ground of known asymptotic results for near- 
polarized and paramagnetic states and not derived by means of general physical arguments. 
In spite of this, also in the intermediate regimes of temperature and magnetic field, this 
formulaprovided results in very good agreement with the exact numerical transfer matrix 
(TM) lO, In] ones for a spin chain and with the exact high-temperature-series (HTS) data 
of Rushbrooke et al. [l2|, for the three-dimensional model. 

Quite recently 0], a study on a class of spin models based on the classical Heisenberg 
Hamiltonian has provided clear evidence of the effectiveness of the old formula for magneti- 
zation suggested in Ref. and of the CDSM itself to describe magnetic properties in a wide 
range of temperatures, in surprising agreement with high resolution simulation predictions. 
On this grounds, the authors were also able to explain some puzzling experimental features, 
confirming again that this formula appears to work surprisingly well in different contexts. 

In this paper we are primarily concerned with this relevant question by accounting for an 
extension of the quantum Callen method to derive, in the context of the classical two-time 
GF-framework, and, on the ground of general arguments, an expression for magnetization 
of the d- dimensional classical Heisenberg model, which is just the classical analogous of 
the famous quantum Callen formula. Remarkably, the formula here derived reproduces 
the asymptotic results obtained within the CSDM in the near-polarized and near-zero 
magnetization regimes. 

To test the most reliability of our formula, here we limit ourselves to calculate relevant 
quantities, as the spontaneous magnetization and critical temperature for d = 3 and other 
ones ioY d = 1, and to compare our predictions with those obtained in Ref. [2|. Noteworthy is 
that, very small deviations are found in the intermediate regime corroborating the accuracy 
of the formula for magnetization suggested many years ago and, in turn, the robustness 
of that found here. 

The paper is organized as follows. First, in Sec. [Ill we introduce the classical spin 
model and the appropriate classical Callen-like two-time Green function with the related 
equation of motion which are the main ingredients for next developments. Besides, we 
introduce, in a unified way, the Tyablikov- and Callen-like decouplings for higher order Green 



3 



functions. In Sec. Illll we extend the quantum Callen approach to derive the general formula 
for magnetization valid for arbitrary dimensionality, temperature and applied magnetic field. 
It is easily obtained overcoming the intrinsic difficulty related to the feature that the classical 
analogous of the quantum spin identities used by Callen ol does not exist. In Sec. IIVI 
self-consistent equations are obtained allowing to determine the magnetization and hence 
other thermodynamic quantities of experimental and numerical interest. Calculations for 
magnetization, transverse correlation length and critical temperature for d > 2 and different 
lattice structure are presented in Sec. |Vl Finally, in Sec. IVl| some conclusions are drawn. 



II. THE MODEL AND CALLEN-LIKE GREEN FUNCTION 

The classical Heisenberg ferromagnet is described by the Hamiltonian 

^ N N 

7/ = -- ^ JijSi -Sj - h'^ Si 

i,j=l i=l 
N N 

= - 2 E '^^^ i^tS- + S^S^) - h J2 (1) 

i,j=l i=l 

Here is the number of sites on a (i- dimensional lattice with unitary spacing, 
{Sj = [Sj, Sj, Sj) ; i = 1, 2, A^} are classical spin vectors with |Sj| = S, = Sj ± iSj, 
Jij = Jji, {Jii = 0) is the spin-spin coupling and h is the applied magnetic field. Of course, 
the identity = (^SjY +Sj' Sj^ = S"^ is valid. For formal simplicity, in the next developments 
we will assume 5=1. This assumption is perfectly legal in the classical context. 

The model ([1]) can be appropriately described in terms of the 2N canonical variables 
(f) = {(f) j} and = {Sj}, where (pj is the angle between the projection of the spin vector 
Sj in the x — y-plane and the x-axis. 

The basic spin Poisson brackets are 

{S!,Sp} = T^StS^,, 

{S+,Sr} = -2zS^6,„ (2) 
{sr,S^} = efS]6.,, (a,/3,7 = x,2/,^), 

where {...,...} denotes a Poisson bracket and e"^ is the Levi-Civita tensor. 
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In strict analogywith the quantum counterpart 9(], we now introduce the classical two- 
time retarded GF [l|, |5| 

G,,{t-t') = e{t-t'){{S+{t-t'),e'''^ST}) 

= {{Sl{t-t'y,e'^'^Sr)). (3) 

Here 6{x) is the usual step function, a is the Callen-like parameter, (...) = 
Z-^ rij lo'' d(pj j\ ...c/S'|e"^^(^'^^^'{'^^'}) stands for the usual statistical average, {A, B} de- 
notes the Poisson bracket of the dynamical variables ^(0,5"^) and B{(f),S^), and X{t) = 
e*^*X where L... = i {Ti, ...} is the Liouville operator. Of course e*-^* acts as a classical time- 
evolution operator which transforms the dynamical variable X = X(0) = X(0(O), 5'^(0)) at 
the initial time t = into the dynamical variable X{t) = X{(f){t), S^{t)) at the time t. 
The equation of motion (EM) for the GF ([3]) is given by (with t = t — t') : 

^^ = Sir){{St,e'^'^Sr}) + {{{srir),n};e'^'^Sr)), (4) 

which, in the frequency-o; Fourier space, becomes: 

ooG,,{u:) = z{{St,e'^'^S-})+z{{{Snr),n};e'^'^S-))^ (5) 

with = {{St{Ty,e'^'^Sr))^ and ((A(r);i?))^ = rfre--((A(r); 5)). 

Now, from basic Poisson brackets ([2]), employing a straightforward algebra yields: 

t{{St,e'''^Sr})=^{a)6,,, (6) 



where 



and 



^(a) = i{{S+,e^'^'Sr}) = -aQ{a) + 2Q'{a) + aQ"{a), (7) 



n(a) = (e'^^0- (8) 
On the other hand, in Eq. (jlj), we easily have also: 

{SI, 'H}=^J2 [S!S,^ - S^S(\ - ihSf. (9) 

h 

Then Eq. ^ becomes 

(w - h)Gij{u) = tp{a)5ij + 
- JiH m{r)S^{Ty, e'^'^S.))^ - ((5+(r)5^(r); e'^^.^^")).] . (10) 



It is worth noting that, for magnetization per spin cr = {S^) we have: 

a = i^(0). (11) 

Up to this stage, no approximations are involved. 

For our purposes, to close Eq. ffTOj) we consider here the following decouplings for the 
higher order-GF's: 



{{SfXr)St{ry,e'^'^S-)) 



a 



(12) 



where A = (TD) and A = 1 (CD) correspond to the Tyablikov- and Callen-like decouplings. 
At this level of approximation, Eq. (fTOl) reduces to the closed equation: 

{u - h)Gij{uj) = ij{a)6ij -crY^h Jih [Ghj{<^) - Gij{uj)] + 

+Af -hh [{SrSt)G,,{u) - {S^St)G^,{uj)\ . (13) 

We now define jol, the k-wave vector Fourier transforms in the first Brillouin zone as: 

^^^■(^) = ]^Y.^''-^"~^'^G^{ujl 

k 

= ^E^^''^""''^^^(k), (14) 

k 

a,(u;) = (5+5r) = l5^e^'^-(--^)C(k). 

k 

By replacing f[T^ in f[T^ . in the (k, a; )- Fourier space the resulting equation for G\^{uj) can 
be immediately solved and we find: 

G,(.) = (15) 



where 
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uj^ = h + aim - ^(k)) + ^ 5^ [J(kO - J(k - k')] C(k'), (16) 

k' 



defines the dispersion relation for the undamped classical oscillations in the systems in terms 
of the Fourier transform J(k) of the spin exchange coupling Jij. 

At this stage, in order to determine G'k(i^), and hence the relevant magnetic properties 
of the classical HM as functions of T and h at a = {] (as for instance a = 1^(0) /2), one must 
obtain an explicit expression of ip{a) in terms of the basic correlation functions related to 
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the original GF. A direct calculation of a constitutes an intrinsic difficulty for classical spin 
models since the classical counterpart of the quantum kinematic rule for the z component 
of the spin does not exist. 

Here, we extend to the classical HM the well known Callen method for calculation of ip{a) 
and hence of magnetization which is the key of the next thermodynamic analysis. This will 
be the main subject of the next section. 

III. CALLEN-LIKE APPROACH FOR MAGNETIZATION 



We ffist note that, from Eq. (fT5|) and the known relation for the spectral density Ak(w), 
corresponding to the Green function G\i{uj) [5], 

Ak(a;) = I [G^{uj + it) - G^{uj - le)] , e ^ 0+, (17) 

we find 

k]^{uj) = 2ti^ {a) 5 {uj -uj]^). (18) 

Then, it is easy to derive the spectral density Ajj(cj) for Gijiu) by using the Fourier repre- 
sentation 

k 

= 27r^(a)^ 5^e^'' (^-'^)5(u; - c^k). (19) 

k 

Now, we can calculate the correlation function {e°'^^SjSf) associated to the initial GF 
G,,{uj) = {{St{T)-e'^''^S-)% 0. 

From the classical spectral theorem |5| 

(B^).rr|^^. (20) 



we get 



/ + 00 
-oo 



duj Ajj(ci;) 
2ti (jJ 



1 . _ gik-(r,-rj) 
k 
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which imphes also that, in Eq. ( !T6|) . 

C(k) = (SjS:,> = ^ = T^. (22) 

In particular, for i = j, we have 

(e'^^»^5r5+) = <|.^(a), (23) 

where the quantity 

does not depend on the Callen-like parameter a. On the other hand, by using the identity 
S^S^ = S^Sl^ = 1 — (S-)^, the correlation function on the left side of Eq. ( 123|1 can be 
independently expressed in terms of the function Q{a) = {e"'^^) as 

{e^''SrSt) = n{a)-n"{a). (25) 

Then, with Eq. ([7]), we finally obtain for Q{a) the equation 

Q"{a) + 2 (^^ + ^'{a) - Q{a) = 0, (26) 

for which the initial condition Q{0) = 1 is valid by definition. This is, of course, insufficient 
to find the physical solution of Eq. (!26l) and one should add a supplementary condition 
to be searched properly. Unfortunately, there is no classical analogue of the operatorial 
identity np^_5 (S^ — p) = which is the key ingredient of the Callen approach for the spin- 
s' quantum HM 9|]. In the following, we will show that, at our level of approximation, the 
additional condition 

n{a) = [ dS'f{S')e''^\ (27) 



which follows formally from the definition of the canonical ensemble average of the dynamical 
variable e"'^^, combined with f2(0) = 1, allows to determine univocally VL{a). 
In view of the structure of the differential equation (|26ll for f2(a), we assume 

f{S^) = giSne''/", (28) 

so that, we can write 

n{a) = J dS' (?(5^)e(*+")'^^ (29) 
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It is worth noting that, in terms of y = l/$ + a {Q{a) <^=^ by introducing 

p{y) = y^iy) Eq. fl26l) reduces to p"(y) — piy) = which can be simply solved by means 
of a combination of hyperbolic sine and cosine functions. So, one easily obtains the general 
solution Q{y) = Asmh.{y)/y + Bcosh.{y)/y where the constants A and B have to be deter- 
mined by using the initial condition = 1 at a = and the average-value nature of 
^{y), as expressed by the additional condition (1271) . Indeed, it is immediate to see that 
the odd solution cosh.{y)/y is incompatible with this last requirement (it should require a 
non-positive weighting function g{S^) introduced in Eq. ( 129|) ). so that it must necessarily be 
B = 0. Besides, from f2(l/$) = 1 one easily obtains A = (1/$)/ sinh(l/$) and and hence, 
as a function of the Callen-like parameter, the physical solution of the original Eq.(26) takes 
the form: 

^ ^ l/<l> + a sinh(l) ^ ^ 

This is the central result of the paper which constitutes the classical analogue of the quantum 

n 

Callen formula [9|. It provides for the magnetization a the expression (which is valid for 
any d, T and h) 



a 



coth ( ^ 1 - $ 



(31) 

where L{x) = cothx — ^ is the Langevin function. Here $ is given by Eq. in terms 
of the classical oscillation spectrum which depends on a itself. Hence, Eq. f l3T]) is a 
self-consistent equation for a where T and h enter the problem through 



IV. SELF-CONSISTENT EQUATIONS 

For next developments, it is convenient to work in terms of the dimensionless quantities 

7k = 7M' = 7^, f = h = j^. (32) 

Then, the basic equations to determine the spectrum of elementary excitations and the 
magnetization become 

= /i + a(l -7k)i?(k) 



with 



(33) 

a = L(l)=coth(i)-$ 

^ ' yiBz(27r)'^(l-7k)a;k'' ^ ^ 
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and 

^-t[ (35) 

Although previous equations are true for generic center-symmetric exchange interactions, 
we focus here on short-range interactions with J (k) = J e^^'^ and hence 7k = ^ ^^^^ 
where z is the coordination number and 5 denotes the nearest-neighbor spin vectors. 
In such a case, it is simple to show that jol 

(7k' - 7k-k') ^ (k') = (1 - 7k) Iv!^ (k') • (36) 
k' k' 

Then, the basic equations become 



'^k = /i + CI (1 - 7k) -R 

Mi) 

with 

d'^k 1 



(37) 



$ = T / , (38) 

JiBZ {2Tif h + a{l--fk)R 

where now R is independent of k and given by the equation 

R=l + \af [ ^ . (39) 

JiBZ {27cy /i + a (1 - 7k) i? 

Of course, this quantity depends on the decoupling used to close the equation of motion for 
the original GF (A = for TD and A = 1 for CD). 

In the limit /i -> 0+ with a > 0, Eqs. (EZD-dM]) reduce to 

Wk = o- (1 - 7k) R 

'fQ(f)\ , (40) 



a = coth 

tq(t] 



with 



TQ (t) 

$ = (41) 



Q (f) = (42) 

and _ 

i?=l + A^Frf(-l). (43) 
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Here 



FA-V. 



d'^k 1 
IBZ (271)" (1 - 7k) ' 



(44) 
(45) 



'iBZ (271)" (1 - 7k) 

from the conventional notation for the so-called "structure sums" Fa (n) = jj J^k (1 ~ 7k)" 
depending only on the lattice structure of the spin model. 

From Eq. (H3l) one immediately obtains the physical solution for R {> 1): 



1 + Vl + 4ATF,(-1) 



(46) 



With this expression, the basic equations for the zero-magnetic field problem assume the 
form 



where 



.(f) 



coth 



1+ ^/l + 4ATFrf(-l) 

ct(T) \ / fQ{f) 



;i-7k) 



TQ(T 

2FA-1) 



<t{T\ 



l + A\TFA-l) 



(47) 



(48) 



Notice that for A = (TD) the quantity Q does not depend on T. 

Previous equations show that the zero-magnetic field problem reduces to solve the single 
self-consistent equation (HOl) for cr ^T^ for dimensionalities for which the integrals (jH]) and 
f H5|) converge and hence long-range order may occur. 

Now, we have all the ingredients to explore the relevant thermodynamic properties of our 
classical model for different values of dimensionality. This will be the subject of the next 
section. 

It is worth emphasizing that previous equations differ from those obtained in Ref. {2I only 
for the expression of magnetization. The use of this formula avoids making the assumptions 
which were not sufficiently justified and tested almost thirty years ago in the context of the 
CSDM. 



V. MAGNETIZATION CALCULATIONS 
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By inspection of Eqs. (I33|) - (l35|) or (137|) -( 139|) . we see that at zero temperature it is $ = 
for any d and hence a = 1 for any value of > (fully polarized ground state). On the other 
hand, for finite temperature and /i = 0, since 1 — 7k — ^k"^ as A; — > (J(0) — J(k) ~ \J iS^) k^ 
as A; —7- 0), the integrals flH|) and f HS]) converge only for d > 2 where long-range order (LRO) 
is expected, consistently with the classical version [13] of the Mermin- Wagner theorem 1^ . 
For d < 2 one has $ = oo at /i = and no LRO may occur at finite temperature, i. e. 
a ^T^ = 0, defining a paramagnetic phase. In this asymptotic scenario, d = 2 assumes the 
role of lower critical dimension for the classical Heisenberg model. 

In the low-temperature regime, close to the polarized state, we have $ ^ 1 ($ — > as 
f ^ 0) so that coth (i) ~ 1 + O (e-^/*) and Eq. dSZD provides a ~ 1 - $ + O (e-^/* ) . So, 
Eqs. (I371)-( l39|) reduce exactly to those explored in detail in Refs. j2-5| 
consider again this regime. 



^ and here we do not 



Besides, at high temperature we have $ 3> 1, which implies a <C 1, and the asymptotic 
results derived in Refs. {2, 5| are reproduced. 

Then, the properties for intermediate values of $ are of main interest for us since, in this 
crossover regime, the expression for cr obtained in the present paper (see Eq. ( 133|) ) differs 
analytically from the corresponding one (in our notations) a = {77^^}^^^ suggested in Ref. 
2] within the framework of the CSDM. This was successfully employed in the recent study 
7| about more complex classical Heisenberg models obtaining accurate results in surprising 
agreement with experiments. 

In the following part of this section, we limit ourselves to consider the cases d = 1 and 
d = 3foT which suitable results exist [2|,Q for a meaningful comparison with our predictions 
based on Eq. (133|) . 

A. The classical ferromagnetic chain 

As mentioned before, the one-dimensional model does not exhibit LRO (a (T) = in zero 
magnetic field) and the quantity of theoretical and experimental interest is the paramagnetic 
susceptibility x = ^^^h^q | • 
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In the asymptotic regimes T ^ 1 and T ^ 1 with $ ^ 1 we find, respectively, 

+ f >1, 



X 



3T 1 ^ 3T ^ 3 



1^ 
3T 



1 2 



(49) 

r< 1, 



:. I2I in the same 
15| with i and 3 



3 

which coincide, as expected, with the corresponding ones derived in Re: 
regime ($ ^ 1). Remarkably, these results differ from the exact ones 
instead of | and | in Eq. (H9ll . 

Similarly, in the regime $ ^ 1, near the polarized state, we find the same results reported 
inRef. Q- 

As already noted at the beginning of this section, only when h ^ and $ is finite, one 
can expect to find differences from our framework and that of Ref. j2j] at the same level of 
decoupling approximation. 

Unfortunately, a comparison of corresponding results can be performed only numerically. 
The reduced magnetization a, derived by Eqs. (P7|) - (|5^ as a function of T and h within the 
CD framework is plotted in Fig. [T] and compared with the corresponding predictions obtained 

1 /2 

jy means of CSDM with cr = { } and with the exact numerical transfer- matrix data 



Remarkably, at finite temperature and higher magnetic fields, our CD results are nearer 
to the exact numerical transfer-matrix data than those found with the CSDM although 
they are very near and essentially coincident at near zero and high temperatures. 

A further signature of the greatest accuracy of the present results at finite values of $ 
arises from numerical calculations for the transverse correlation length defined by 



Ci= — lim 



(50) 



In general, using the expression of {S^S_y.'^ in terms of a, one easily finds within our 
approximation and in d dimensions 

which should provide information also about the paramagnetic susceptibility X = | decreas- 
ing to zero the magnetic field. 

Numerical CD results of Eq. ( 15T]) for d = 1 with z = 2 are plotted in Fig. 2 and compared 
with other data as in Fig. [TJ 
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FIG. 1: Plots of magnetization a as a function of T for different values of the reduced magnetic 
field h. The solid and dashed lines represent the Callen-like decoupling (CD) results and those 
3], respectively. The dots denote the exact transfer- matrix data. 



derived in Ref. 



Notice that our results appear again more accurate at intermediate values of $ reproduc- 
ing more effectively the exact numerical transfer-matrix results. 



B. d> 2 and three-dimensional ferromagnetic models 

For all d > 2, LRO is possible and the zero-magnetic field equations fH6|) - P8|) allow to 
find the reduced critical temperature and the spontaneous magnetization a {t^ in the 
full range < f < T;. 

First, we note that, for a — t- 0^ (with T — t- and to be determined), we have 
$ = — — > oo. Thus, using the expansion cothx — \ = Yl^=i {2n)" where i?2n are 
the Bernoulli numbers, we can write for cr ^ 1 the equation 

1 / a 



a 



3 ^ TQ (t) / \TQ (t 




(52) 
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FIG. 2: Plot of the transverse correlation length ^± as a function of reduced temperature T for 
different values of the reduced magnetic field h. The solid and dashed lines stand for the present 
CD-predictions and the CSDM ones of Ref. P], respectively. The dots denote the exact transfer- 
matrix data. 



Then, we get 

and Tc is determined by the equation 



1-3T,Q TJ =0. 



(53) 



(54) 



So, from Eq. (|53|) . we have 



1/2 



(55) 



which defines the mean-field approximation (MFA) critical exponent /5 = |. 

The explicit expression of Tc for d > 2 can be immediately determined form Eqs. (jlS 
and (15^ providing 

, AF,(-1) 



3F,(-1) 



(56) 
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TABLE I: Numerical values of the reduced critical temperature Tc of the three-dimensional classical 
Heisenberg model for different lattice structures. HTS stands for the exact high-temperature-series 
results. 





MFA 


CD 


TD 


CSDM 


HTS 


sc{z = 


6) 


0.333 


0.245 


0.220 


0.245 


0.241 


bcc(2; = 


= 8) 


0.333 


0.262 


0.239 


0.262 


0.257 


fcc(z = 


12) 


0.333 


0.269 


0.248 


0.269 


0.265 



In particular, we get 

fico, _ fiTO, L ^ > f<™> (58) 

similarly to the Callen result for the quantum Heisenberg model j^. 

The numerical values of for the three-dimensional Heisenberg ferromagnet, as obtained 
by means of different methods, are reported in Table I for simple cubic (sc) - body centered 
cubic (bcc) - and face centered cubic (fee) spin lattices. 

Notice that our CD-results, as the identical CSDM ones, are very close to the exact HTS 
results of Ref. [12 1. 

As concerning the spontaneous magnetization a ^T^ as a solution of the self-consistent 
equation f H7|) . one expects that the CD results will differ from the CSDM ^ ones in a 
temperature range sufficiently far from the asymptotic regimes near zero ($ ^ 1) and 
critical ($ <^ 1) temperatures. 

Eq. ( H71) can be solved only numerically for any temperature in the interval < T < 
Tc. The data for a ^T^ of the {d = 3)-classical Heisenberg ferromagnet for different lattice 
structures are plotted in Fig. 3 and compared with the CSDM and MFA results. 

It is worth emphasizing again that, as expected, differences between our results and the 
CSDM ones occur for finite temperature in the range 0.05 ^ T < 0.25. 
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VI. CONCLUSIONS 



The main purpose of the present paper was the derivation of a new formula for the 
magnetization of the c?- dimensional classical Heisenberg ferromagnetic model. This is a 
central problem within the two-time GF framework and the strictly related spectral density 
method in classical statistical physics. The main difficulty indeed is crucially due to the 
absence of the classical counterpart of the exact kinematic rule for the 2;-components of the 
quantum spin vectors. By a suitable extension of the Callen method developed almost five 
decades ago for a spin-S quantum Heisenberg ferromagnet js], this new formula has been 
obtained analytically overcoming the previously mentioned difficulty. Surprisingly, although 
it is formally different from that suggested within the CSDM, both the frameworks provide 
results in good agreement with the exact numerical transfer-matrix data available for the 
classical Heisenberg ferromagnetic chain and with the exact HTS ones for three- dimensional 
case. However, the present formula provides more accurate results in an appreciable range 
of temperatures and magnetic field where the differences between the two formulas become 
evident and suitably interpolates the low- and high- temperature regimes. This emerges 
from a careful comparison between our predictions and those numerically accessible for 
the magnetization and other directly related static quantities as the transverse correlation 




0.1 0.2 0.3 0.4 




0.1 0.2 0.3 0.4 



FIG. 3: Plots of magnetization o" as a function of the reduced temperature T for different three- 
dimensional lattices: (a) simple cubic (sc); (b) body centered cubic (bcc); (c) face centered cubic 
(fee). Comparison is made between present TD and CD calculations and the CSDM and MFA 
corresponding ones. 
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length and the paramagnetic susceptibihty. In any case both formulas work surprisingly well. 
Of course, other experimentally relevant thermodynamic quantities, such as the internal 
and free energies and the specific heat at fixed magnetic field, can be found by using the 
general classical two-time GF formalism introduced by Bogoliubov Jr and Sadovnikov l| 
and further developed in successive papers j2"5|, in strict analogy with the better known 
quantum counterpart, but these calculations are beyond the purposes of the present paper. 

In conclusion, we believe that our new formula for magnetization may be usefully em- 
ployed for further theoretical and numerical calculations to explore the properties of other 
Heisenberg-type classical spin systems at different dimensionalities. 
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